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REGULAR SINGULAR POINTS OF A SYSTEM OF HOMO- 
GENEOUS LINEAR DIFFERENTIAL EQUATIONS OF 
THE FIRST ORDER* 

By Otto Dunkel. 

Presented by Maxime B8cher, May 14, 1902. Received July 10, 1902. 

We will consider in the present paper a system of n differential equa- 
tions of the form : 

* 

in which the /i tl /a are constants, and the a,,/s are functions, not neces- 
sarily analytic, of the real independent variable x, continuous in the 
interval : 

< x < b.f 

We shall require that \a t j\ be integrable up to the point x = 0. 
For the development of certain sets of solutions we shall make the 
additional requirement that even after certain powers of log x have been 
multiplied into | a {J | the resulting function shall be integrable up to the 
point x = ; this requirement will be stated more explicitly later. J 

The point x = satisfying these conditions may be called a regular 
singular point of the system of equations (1) in conformity with the use 
of that term by Professor Bocher in the study of linear differential 

* This paper was accepted in June, 1902, by the Faculty of Arts and Sciences 
of Harvard University in fulfilment of the requirement of a thesis for the degree 
of Doctor of Philosophy. 

t The requirement that the functions a,, should be continuous in < x < b is 
made only for the sake of simplicity. We might allow them to have a finite 
number of discontinuities in b of such a kind that each function | a i . | can be 
integrated throughout the interval ; and all the following work would hold with 
very little change. 

% Cf. p. 367. 
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equations of the second order.* That this terminology is a legitimate 
extension of that commonly used when the coefficients of the system of 
differential equations are analytic functions of a complex variable, will be 
evident if the results of the present paper are compared with the thesis 
by Sauvage : Theorie generale des systemes d' equations differentielhs 
lineaires et komogenes.f 

Our object is to investigate the nature of the solutions of (1 ) in the 
neighborhood of the regular singular point x = 0; and for this purpose 
we shall first reduce the system of equations to a canonical form by 
means of a linear transformation with constant coefficients of the de- 
pendent variable. We shall then apply the method of successive 
approximations to develop about the point x = a system of n linearly 
independent solutions of the canonical system. By means of the linear 
transformation we shall return to n linearly independent solutions of the 
original system ; and finally an application will be made to the case of 
the single homogeneous linear differential equation of the «th order. 



§1. 

A Special System op Equations: Its Reduction to a 
Canonical Form, and Solution, t 

Let us first examine the special case of (1) in which the coefficients 
a, j are all zero. In this case we have the system of differential 
equations : 

(3) 1=2^ (-1,V..„), 

where the /a,-,/s are constants. 

A solution of this system may be obtained in the following way. 
Substitute 

y { = C,x r C ( = constant 



* Cf. Trans. Am. Math. Soc, Vol. I. Jan. 1900, p. 41. The results of this paper 
are included as a special case in those we now give. Cf. § 7. 

t Paris, 1895. Reprinted from the Annales de la Faculte des Sciences de 
Toulouse, Vols. VIII. and IX. 

I The results of this section are not new, being on the one hand only slightly 
modified forms of Weierstrass's results (cf. the foot-note on p. 345), and on the 
other hand special cases of the results obtained by Sauvage (cf. the last foot-note.) 
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in the equations (3), and then determine the constants r and C { so that 
the equations are satisfied. We obtain in this way the following system 
of n homogeneous linear equations for the C's : 

(4) ft 1 C 1 + ... + ( ft ;,-f)C+... + ft„(?.= (i = l,2,. ..n> 

The necessary and sufficient condition that (4) may be satisfied by a 
set of O'a not all zero is that the determinant: 



(5) 



A(r). 



Mi, i 



M«,i 



F-i, 



/*«, n 



shall vanish. This determinant equated to zero gives an equation of the 
nth degree in r, which is called the characteristic equation of (3) ; * the 
determinant itself we may call the characteristic determinant of (3). 

If the characteristic equation has n distinct roots r u r 2 , . . . r n , we can 
determine n linearly independent solutions of (3) of the form : 



(6) 



9u =<><.,# 



i — 1, 2, . . . n 
i = 1, 2, . . . n 



If however there is a multiple root, there will be, in general, solutions 
involving powers of log x. To determine these solutions, we must 
examine the minors of the characteristic determinant (5), and ascertain 
if this multiple root is also a root of all the first minors, second minors, 
etc. For the further study of this case, it will be useful to introduce the 
conception of the elementary divisors of A(r). 

Suppose r' is a root of the determinant A (r) such that all the pth minors 
of A(r) are divisible by (r — r 1 ")*, but no higher power of r — / divides 
them all. In the same way (r — r ,s ) p + 1 shall be the highest power of 
r — r' dividing all the (p + l)th minors. Then the expression : 



(r - r'f 



= 1,-1 



/+i 



* Sauvage : I. c. p. 80. 
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is called by Weierstrass an elementary divisor of the determinant A (r).* 
It will be convenient to employ a different notation from that used in 
the definition of an elementary divisor. An elementary divisor of A (r) 
will be written : 

<r-r t y; 

and it is to be noticed that several r's with different subscripts, may be 
equal, as will be the case when a multiple root furnishes several elemen- 
tary divisors. We shall always have : 

k=m k=m 

(7) A(r)=Yl( ri -ry* ^e t = n. 

It can be shown (cf. the next foot-note) that a necessary and sufficient 
condition that a pair of systems of differential equations : 

, j—n # j j=n 

dx ~ f* x yj ' dx ~ ** x ' 

(t = 1, 2, . . . n) 

can be transformed the one into the other by means of a transformation : 

j=n 

y, = 2^,. y zj (i = 1, 2, . . . n) , 

in which the -4's are constants whose determinant is not zero, is that the 
characteristic determinants of the two systems have the same elementary 
divisors. This theorem enables us to simplify the solution of the system 
(3) ; for we can write down a second system of differential equations 
having the same elementary divisors as (3), as follows : 

(8) &-^* w+ ***' {l=l,2,...e t )' 
where z*„ = 0. 



* Cf. Muth : Theorie und Anwendung der Elementartheiler, p. 2. 
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(9) A(r) = 



r x — r 
1 r x — r 



1 r x — r 



r 2 — r 
1 r 2 — r 



1 r 2 — r 



r — r 

1 r m — r 

1 r 



and it will be easily seen that it has the elementary divisors (r — r-^f 1 , 
(r — r 2 ) e \ ... (r — r,^f m . Then, by the theorem above referred to, 
there exists a set of n 2 constants A ikl) whose determinant is not zero, 
such that: 



(10) 



k=l 1=1 



(» = 1, 2, . . . n). 



The system of differential equations (8) we may speak of as the 
canonical system ; and now it is easily seen that : 

The canonical system of equations admits e K solutions, corresponding to 
the elementary divisor (r — 7V) e "> of the following simple form .- 



(11) z Kt i=0 



«K,l = 



Q-xy. 



x r " (log x) 



l—k 



I < A. 
A. < I < e K 



(A = 1, 2, . . . e K ) 



and the n solutions obtained by giving k the values 1,2, ... m are linearly 
independent. 

* Cf. Weierstrass, Werke, Vol. II. pp. 75, 76. The case considered by Weier- 
strass is very easily reduced to the one we are considering by the change of inde- 
pendent variable t = log x. This reduction of Weierstrass is also given in Muth's 
Eleinentartheiler, pp. 195, 198. On page 198 are a number of references to the 
use of the theory of elementary divisors in the study of differential equations. 
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The determinant of these n solutions can be written out in such a way 
that the elements above the principal diagonal are all zero. Thus the 
value of the determinant is : 

K—m 

k:=1 

X , 

and the n solutions are therefore linearly independent. 

On account of the relation (10), each solution (11) of the system (8) 
will determine a solution of the original system (3). Accordingly we 
have n solutions of (3) which are linearly independent, for their de- 
terminant at any point is equal to the determinant of the A's in (10) 
multiplied into the determinant of the solutions (11) for the same point, 
and neither of these determinants is zero. 

Suppose now we consider any multiple root of the characteristic de- 
terminant (5) ; for simplicity let us take r u and suppose that 

(12) r x - r 2 = . . . = r k , 
so that 

(13) e, + e 2 + . . . + «, 

is the multiplicity of the root r 1# Then from (10) and (11) we see that, 
corresponding to this root, there are h solutions of (3) not involving log x 
and linearly independent : 

(») ».- = A,,*" (.':j;J ;;;;). 

Therefore the constants A i K e are linearly independent solutions of the 
equations (4) when r = r x , as we readily see by putting the values of 
y/'"" in (3). Now the equations (4), in this case, have only k linearly 
independent solutions, since there are only h elementary divisors corre- 
sponding to r x , and therefore at least one &th minor in (5) is not zero 
when r = r x . If, then, we have any other solution of (3) of the form : 

y,= C,x r > (.' = 1, 2, . . . «), 

the constants C £ must be linear combinations of the k sets of constants 
-4,- , „ • We have then the following result : 
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The system of equations (3) admits n linearly independent solutions, 
such that, corresponding to each elementary divisor (r — »v) e " of the 
characteristic determinant, there are e K solutions : 

(") ^ = * r *lVx)! *-« (,og * y-A (x = 1 1 : : : 1) • 

If r K is a multiple root of the characteristic equation which furnishes s 
elementary divisors with the exponents e K , e,^, . . . e, s _ 1 , then the 
constants : 

A i,K,e K > A i,K+l,e K+1 > ■ • ■ J i,K+s-l, Vf s -l' (* = 1, 2, . . . n) 

are s linearly independent solutions of the equations (4) when r = r K . 



§ 2. 

Solution of the Canonical System in the General Case 
bt Successive Approximations. 

We shall now return to the system of equations (1) ; and here again 
we shall make use of the linear transformation (10) to reduce the system 
to the canonical form : 

< 16 > £**' = x •*« + S**< +22 *Z*u * 

ax x x i=1 ^._ 1 

(issl, 2, ...«) (Z=l, 2, . ..«. 

The coefficients 6^ are linear functions with constant coefficients 
of the coefficients a t/ in (1). 

We shall now make use of the method of successive approximations to 
develop solutions of (16) about the point x — 0. It will be convenient 
to write the equations (16) in the form : 

The first approximation will be indicated by a third subscript 0, and is 
obtained as a solution of the system of equations resulting from (17) by 
making the right side zero : 

* This reduction is used by Sauvage in the case of a system of equations with 
analytic coefficients. L. c, pp. 89, 90. 
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na\ d \ r k 

( 18 ) dx Zi ''- a ~ x**-'- 1 ' ~ x Z *-''° = 

We have seen how to solve this system. Any one of the n solutions 
that we have obtained may be used as the first approximation. This 
approximation having been chosen, we insert it in the right side of (17) 
and obtain the following relations for the first correction : 

d 1 r k ^ v;'« 

ax x x i _ x j _ 1 

The right side is now a known function of x; and we have, conse- 
quently, a system of non-homogeneous linear differential equations to 
solve for s^ t< 1 . Having determined this first correction, it is inserted in 
the right side of (17), and the resulting equations are solved for the 
second correction. This process is repeated again and again, the relation 
connecting the <?th and the (q + l)th correction being : 

, , i=m r—e, 

,ia\ d 1 r » - X ^ V' y 

V > jZ Z *.l<!+ 1 ~ ~**.'-l.S+l ~~ ~Z Z ^',l+>- — £t Z* °k,l Z i.j,!' 

ax x j. i _ 1 j=l 

Each equation (19) may be written: 



x *^z( x r * z *.i,v+i)-~ z *.'-i.!+i= 2 2 h Xt z <.J.<> 
"* x t=i j=i 

whence : 



( 2 °) 

**, i. ,+l = * r »\ J x ~ ~^ k *k. /-I, rt-i dx + J x-** 2 2 Ki z 'J. i dx \- 



k,l "£,l 



Now writing out the value of Zj.,,_i j?+ i in the same way, and substitut- 
ing it in the first integral of (20), we have : 

xx . x 

- dx / x- r *^ 2 b li-x z <j,* dx + / a5-r *2 2 *i / *«»«*» • 
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Now substitute in this result the value of ^,/_ 2 , ? +i ; in the result thus 
obtained the value of z kJ _ Si g+1 , etc. After a certain number of substitu- 
tions we have : 

X X 

(21) z, i! ,^ = x r *\j\dx...jx- l - r *z t . Lig + l dx 

C k,l C i,L+l 

+ ^j-dx...\ x+*2 ^ b'^z^dx , 

where L < I. 

When L = we have : 

X X 

(22) z*,,,^! =x r *2 ml I -dx . . . / x~ r » 2 2 %n-i-t*tj.t dx - 

c k, I c t, i+l-l 

The lower limits e tt , will be determined later to satisfy several con- 
ditions. 

We may choose at pleasure any one of the m elementary divisors, say 
(r — r K ) e,c , and then select any one of the corresponding e K solutions of 
(18) for the first approximation. We shall take then for the first ap- 
proximation z^/o> the values given in (11) for a particular A; the inte- 
gers k and A will remain fixed for the solution we are now developing. 

For the development of the solutions corresponding to (r — r^) e ", we 
shall make the following further assumption as to the coefficients b'ifp 
Let us examine all the exponents of the elementary divisors (r — r^) e *, 
which are such that Br k = Rr K , where Rr k means "real part of r^.," 
and pick out one exponent, say e K , that is as great as any one in this 
special set of exponents ; i.e., 

(23) e K > e k , where Br k = Rr K . 
The assumption is that the integrals : 

I 



/icin».f-* (i=t?. 



converge. If in particular r K is a simple root, and no multiple root has 
the same real part as it, then e K — 1 = 0, and this further restriction 
drops out. Or it might happen that r K is a multiple root, but that all the 



I. 


Rr k > Rr K 


c k,i = c 




II. 


Rr k = Rr K 


{e M = 


I 51 Z = e^- — e K + A 
Z>Z 


HI. 


Rr k < Rr K 


c 4 ,/ = 
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exponents of the set (23) are unity, and in this case the restriction drops 
out also. 

The lower limits of integration in (22) will be determined as follows : 



(24) 



where c is a constant not zero, which will be determined more closely 
later (cf. p. 358). It will be proved in § 3, that even in the cases in 
which the lower limit is zero, the integrals converge. 

When all these conditions have been satisfied, we build from the first 
approximation and the successive corrections the n infinite series : 

>? /4 = 1, 2, . . . m\ 

(25) * i=a ^**' Um,-..J 

which will be proved in § 4 to converge and to form a solution of the 
system of equations (16). 

It will be convenient to consider in place of the functions z k lq certain 
new functions <f> k ,i, q , which will be defined by the following formulae: 



x r < </>. 



'*, l, 7 



when 


Rr k -^ Rr K , 


or 


Rr M = Rr K , h -if. k and I < L, 


or 


h = k and I < A, 



(26) g^ t = x r * (log x) l - L 4>^ g Rr t =Rr„ k + k and L < I < e t , 



*«,«,, = ^0°g*) <t>«,i, q X<l<e K . 

For the case of q = 0, i. e., the first approximation, the <£'s are certain 
constants (cf. (11)); for all other values of q the formulae (26) define 
them as continuous functions of x so long as x is not zero ; and we shall 
see later (Proof of Convergence, p. 358 et seq.) that each one approaches 
zero when x approaches zero. We shall therefore define each <£*,/,,, 
when q > 0, as zero for x = ; and with this definition they will be 
continuous functions of x in the whole interval Ob. 

These <£'s can be computed from the following recurrent formulae, 
which are easily obtained from (21), (22), (24), and (26) : 
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I. Rr*> Rr K . 



(27) 4W! = a*"^ - " ^ f-^- 



x 

c 

X 



-jxj *-^> j 2 ^ ^ (log x)k/ dXi 

CO *= 1 J =1 

t integrations 

where h tJ is a positive integer whose value it will not be necessary to 
write out, noting, however, that : , 

(28) K,<e K -\. 

II. Rr^-Rr,,. 

a) k 4= k. I < L . 

(29) ** Ar+ i = af*-"' X f-<&- 




a; x 



t=m ,;=:<•,. 



o ,=1 •'— 1 

£ integrations 

J) £ rf 1 k, L < I <, e k . 

m **«* " s5^[/>- ' -is*/: ""'""■' *«■« * 

c c c 

(£ — X) integrations 



tf integrations 

c) k = k, I < \. 



( 31 ) *«,/,rt-i=2 J x^'" 



■h dx 2 2 £«-«*«» Gog *)**'<&. 



o 6 tel •>• 

if integrations 
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d) k=z K , \<l< L. 

t=l x 

(32) fa, t.i+1 = — -j^ ^ \ -dx ... 

(log*)' tt-J * 

* *i=m j=e,. 

•■•Jx dx J 2 2 *«^n-«*w f Gog*)*^«fe- 

o o *— 1 J— 1 

I integrations 

e) k = k, L < I < e K . 

< 33) i r h h h 

^■ l '^ = ilogx^lJ x dx -'-J x dx J x**.^Qogxy* e *dx 

C C C 

(l~L) integrations 
t integrations 

III. JRr i < Er K . 

Here we have the formula (29) again. 
In § 4 we shall consider the n series : 

(34) **,= *!*«, ( k r]'i'--- m \, 

£o \l=l,2,...e l J 

which are such that if we multiply each hy its proper factor x r * (loga:)**.', 
where h kt is given in (26), we obtain the n series (25). We thus 
reduce the proof of convergence of (25) to the question of the con- 
vergence of (34) , and this last question will be settled by reference to 
certain formulae to be established in the next section. 



§3. 

Lemmas Concerning Multiple Integrals. 

We now prove a number of lemmas, which will be useful in the proof 
of convergence, and which also verify the statement that we have made 
that the integrals of the last section, in which the lower limit is zero, 
converge. 

Lemma I. If b is a Junction of x, continuous in the interval 
< x < c, and such that | b | |.log x |' _I , (t an integer > 1) is in- 



DUNKEL. — LINEAR DIFFERENTIAL EQUATIONS. 353 

tegrable up to x = 0, then b can be integrated t times from x = as 
follows : 

X XX 

(35) f {x) = f-dx...f~dx Cbdx. 





t integrations 



Lemma II. If the conditions of Lemma I. hold, then : 

X XX X 

(36) \f t (x)\ = \f~dx...f~dxfbdx\< f\b\ \logx\'-'dx. 
'o o o ' 6 

t integrations 

(0<x < 1). 

We will prove these lemmas by mathematical induction. They are 
true when t = 1. Let us assume that they are true for a particular 
value of t, say t = t t . 

Let X be any particular value of x in the interval < x < 1, and 
choose e at pleasure such that < e < X. Then we have : 

X X X x 

\C \ f h (x) dx I <J~\ f h (x) \dx<J\ dxj\ b\\\ogx I*'" 1 dx 

c e « 

X < 

= logX C\ b 1 1 log x \ k ~ l dx - log € fl&Hloga;!' 1 - 1 ^ 

x 
— floga:|6| \\ogx\ tl ~ x dx 

e 
« X 

= I log e | C\ b | | log x I' 1-1 dx + C\ b 1 1 log x I' 1 dx 

« 

X 

-|logX| flbWlogxf^dx 
o 

« X 

^ |l°g <\f\l\ llogzf 1-1 ^ +J\ b\ I log x^dx 

5 « 

< x x 

< f|S| | log a; I' 1 dar + f\ b \ | log x \ h dx = f\ b | | log x f> dx. 



o 

vol. xxxvin. — 23 
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Therefore, when « approaches zero, the integral on the left side of all 
these inequalities converges, and we have : 

x i x 

lA+i (*) I = IJ'Ia (*) d * I ^/l * 1 1 lo s * I' 1 & - 

o o 

Thus the two lemmas are proved when x < 1 ; and it is easily seen how 
to conclude the proof of I. in case c > 1. 

Lemma III. If b is a continuous function of x in the interval 
< x ^ c, ararf tfc absolute value is integrable up to x '= 0, and if: 

X XX 

F,(x) — I -dx . . . j -dx I x r bdx 


t integrations 

where r is real and greater than zero, then : 

x 
(37) | *-*•,(*) | ^J"| 6 1 dr. 

6 

When t == 1 we have : 

« X X X 

I «-* F x (x) I = af- 1 C 3? b dx < x~r Caf \b\dx < f\b\dx, 
o oo 

and so in this case III. is true. Assume that it is true for t = ^ ; then 
it is also true for t = t x + 1. For: 



■F^(x)\ = x-\p-F h (x)dx <x-f-\F h (x)\dz 
6 o 

x r X 

^^Jx\7^j [h{dx \ dx 



* fi.h* dx ] [/' h • *] = 7,/ 1 6 ! *• 



Therefore III. is true for all values of <. 



DUNKEL. — LINEAR DIFFERENTIAL EQUATIONS. 355 

Lemma IV. Ifb is a continuous function ofx in the interval < x <^ c, 
and if: 

X XX 

g, (x) = I - dx . . . I -dx I x~ r b dx 

c c c 

t integrations 

where r is real and greater than zero, then : 

X 

(38) \^ gi (x)\<~J\b\\dx\ (0 <X ^c). 

c 
When t = 1 we have : 

X X 

I offfi 0) I ^ *? fx- r \b\ \dx\ ^ f) b 1 1 dx\ . 

c c 

Assume that IV. is true for t = t x . Then : 

* XX 

I *- 9 h+x (x) \<»ff-\g ti (x)\\dx\< ^if*- 1 - 1" f\ b | [ dx | J \dx\ 

^^[/ af ^ , * sl ][/ ,4ll ' fal ] = ^[ 1 -(0']/ ,i,ldirl 

c c c 

X 

^jj\b\\dx\. 

c 

Therefore IV. is true for t = t t + 1, and the lemma is proved. 

Lemma V. If b, r, and g, are defined as in IV. and it is further as- 
sumed that the absolute value of b is integrable up to x = 0, then : 

(39) limit x r g t (x) = 0. 

o=0 

To prove this, let us choose a constant s such that < s < r. Then : 

X XX 

\ x "9> O) I = x s \x r - s j -dx . . -\\dx fx^-*x- s bdx\ 

C CO ' x 

c 

as we see from IV. by replacing r by r — s and b by x~ s b. 
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If e is a positive number chosen arbitrarily small, we can choose 17 so 
near that : , 

f \b\\dx\<\, 
and then X ^ 17 so that : 



£ 

2 
Then 



X s fx- s \b\\dx\< 

C 

X s \x~ s \b\ \dx\ = af I ar*\b\ \dx\ + a? \ x~*\b\ \dx\ 

C C 1) 

17 X 

<,X' Cx~*\b\ \dx\+f\b\ \dx\ (0<x<X^ v ) 

c ' 

± + f\b\\dx\< t . 



< 



Therefore : 

£ 



and V. is proved. 



\ a rg t {x)\< (r _ sy _ i (0<x<X) 



Lemma VI. If fi is a function of x continuous in the interval 
< x <i c, and 

(40) limit B = ; 

and if: 

X X 

G/ K x)=j\dx...J l -Bdx 

c c 

t integrations 

then : 

(41) So^^- 80 - 

When t = it is obvious that (41) is true, for then : 
limit G (x) = limit B = 0. 

Assume now that (41) is true when t = t x ; then it will also be true 
f or t — t x + 1 . For if £ is a positive number chosen arbitrarily small, 
we can choose t] so near that: 
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357 



| (log x)h 



G k (*) < 5 & + 1) (0 < X < r, < 1) 



X f\ X 

G h+1 (x) = p-G k (x) dx = J l -G k (x) dx + fi G h (x) dx 

C C If 

X X X 

iJi G h (x) dx\<f l -\ G tl (x) | \dx\ < I& + l)f^^\dx\ 

» * * (0 < x < tj) 

e | (loga;y.+l— (logr?)^ 1 ! 

= 2 ( ' i + 1) | i^n I 

= iiiogx|. +1 [i-(^y i+l ]<ino g zi^. 

Now choose X in the interval < x < -q so that : 



Then: 
1 



(log a:) 



1 Jl Gtl (x)dx 



(0<x< X). 



(log a3)*i+ a 



<?(,+! (*) 



< 



(log a:)': 



l^JI^^^I+l^-l^JI^^)^ 



and : 



(0 < a; < X), 



limit - 



1 



G tl +i(x) = 0. 



x=o (log a;) to- 1 
Therefore VI. is true for all values of t. 

Lemma VII. If (i is a continuous function of x in the interval 
< x ^ c < 1, o«rf smcA that it is not greater in absolute value than the 
constant N, and if G, is defined as in VI, then : 

X X 

(42) | G t (x)\ = I f-dx . . . f-/3dx\ < .ST | log a; |'. 

1 1/ X if X i 

c t integrations 
r • x x x x 

\G l (x)\<j\\dx\...j\N\dx\<NJ l -\dx\...jl\dx\ 

c c c c 

<^[/l|^|]'<^|log,r[l-g|J<^|log, r . 



For: 
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§4. 

Proof of Convergence. 

The convergence of the series (34) will now be proved ; and in this 
proof it will also be shown that : 

(43) limit* =0 9 >0. 

xz=Q ' * v 

We shall use in the proof the following functions. A function of x, 
B, is chosen having the same properties of integrability as b\' t , and such 
that: 

w b> 4%m- (?:"::::)■ 

Next consider all the differences Rr k — Rr K which are not zero, and 
choose a positive constant d such that : 

< d< \Rr k — Rr K \. 

Then (J shall be a positive constant such that : 

fc=» -I 

(45) C>n, G^^-^. 

t—i a 
We also define : 

x 

(46) M(x) = CoB^ogxf^dx. 

o 
The point c is chosen so that : 

(47) < c < 6, c < - (where loge = 1), M (c) < 1. 

This is the final determination of c to which we have referred on page 
350; and this point c will, from now on, mark one end of the interval 
for x. Instead of M(c), we shall, for the sake of brevity, write 
simply M. 

The convergence of the series (34) will be proved by showing by 
mathematical induction that the following inequalities hold for all 
values of q : 



* We might, for instance, take B — 2 2 2 2 | b*^ | 



Tcz=.m fee. i=zm j=e t 

,2 2 2 2 | 
fc=i i—i i=i j=i 
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(48) 



i*-is* (f=iS:::;> 



For q = they are obviously true. Assuming them true for a special 
value of q, q — ft, we will consider the cases outlined in (24) in turn. 

I. From (27) we have : 

(49) l* M)?1+1 l 1 

x R(r k -r K ) ^ fl\dx\...f±\dz\f X-^-"^ BM* | log X f «~ A | dx \ 

— c c c 

t integrations 

after replacing 

IW> i^*K 2 21^1 

respectively by the greater values 

M 9 ', \]ogx\ e ~\ B 
(Cf. (28), (30), (32)). 

t-=i -i * 

l*M, 9l+ ll ^ ^^[^Z^J'/^llog^r^l^KLenimalV.) 

C 
X z 

< M gi C fB\\ogx\ e «- x \dx\ < M qi f OB\log x ^ \ dx | 

c 

< M Ql M{x) £ M 9l+1 . 

From Lemma V. it follows that the limit of the right side of (49) is 
zero, when x approaches zero ; for we have assumed that B | log x \ 6sr ~ 
is integrable up to x = (cf. page 358), and therefore B | log x | e " — x 
must also be integrable up to x = 0. So (43) is verified for q = q x + 1. 

II. From (29) we have for case (a) : 

+ ; X XX 

(50) |^,,, ?I+1 |< 2 j\ dx - .J\dxJBM^\\ogx\ e ^dx 

*- 1 

t integrations 

t=i % 

■£M qi ^ I B | log x | e «-*+'-i dx (Lemma II.) 
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since the greatest value of e K — X + t — 1 is 

e K — X + L — 1 = e, - A + (e K — e K + X) - 1 = e K — 1. 

X 

( 51 ) l**,i l9l +i[ = M«LjB\\ogxf^dx = Jf^UX*) 

o 

< ^''.^(a:)^ -¥ <7l+1 . 

It is obvious that the above result holds also for II. (c) (31). From 
(32) we shall obtain in the same way for II. (d) : 

X 

(52) | log x I 1 -* | ^ iUi+1 1 < J/" X J*i? | log z l^ 1 ^ = M«> £ M(x) 

o 

I log * I O 

In all three sub-cases (a), (c), and (c?), it is easily seen that (43) is 
true for q = y t + 1- We have now left of case II. the sub-cases (b) 
and (e). 

From (30) we have for case (b) : 

(i — L) integrations 
/— 7 T x x x 

+ 2 J^i...J , ^i«fci/^^ s, iiog*i e "- A i«fei]. 

<=1 c c c 

2 integrations 

In the first part of the bracket we have replaced | <f> k L , a | by 
M qi - 1 M(x), using inequality (51). 

We will consider the two parts of (53) separately. 
(54) x xx 

jl _^/I|^|.../I|^|/I^^(,)|^|<i/^^) 

C CO 

{I — L) integrations 

= i/ 9l+1 - (Lemma VII.) 

G 



DUNKEL. — LINEAR DIFFERENTIAL EQUATIONS. 361 

(-) (1 ^p Yfl^l ■ ■ .flwfBM*]***^]* 

c c c 

t integrations 

< 

t- 



? ' l0g X r [ <\ogx\'- L I B ' l0g * ^ ^ (Lemma VIL) 

< jf ft <iz^ 5> f # I l g x |^-i «& < 1/ft+i (till . 
log a; J ' 8 ' - C 



o 
Therefore : 



I <A I < if" 1 " 1 " 1 F~- + / ~ Z 1 - Jf 9l+1 - < J/ft+ 1 

The limit of the left side of (54) for x = is zero by Lemma VI., 
while the same thing is true of the left side of (55) from the inequali- 
ties. Therefore (43) is true in this case for q = q x + 1. 
From (33) we have for case (e) : 

(56) 

£ X X 

I^^IS nS5 ^[/i|*|.../i|&|/iir^(-)|ft| 



G G C 

(I — L) integrations 



t—l—L 



+ 2 fl\te\...fl\dx\fBM«i\\ogx\**-*\dx\]. 

c c c 

t integrations 

In the first part of the bracket we have used the inequality obtained 

from (52) : 

\\ogxf^\^ Lqi+1 \^M^^M{ X ). 

Now the only difference between the inequalities (53) and (56) is in 
the power of | log x | outside the brackets ; and, since ] log x f~ K >. 
| log x \'~ L , all the results that we have obtained from (53) will follow 
also from (56). Then for all the sub-cases of II., (43) and (48) are 
true for q = q 1 + 1. 
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III. Here we have : 



< x-*^-*') ^ f-dx...f\dx /V (r "- r *> B M q >\ log x \ e *-* dx 

*- x 

t integrations 

<M qi 5 1^7^ n ' fB\logx\ e -- x dx (Lemma III.) 

X 

<M qi C 



CB\\og,x\ e K-*dx< M 9l+1 . 
o 

It is obvious that (43) is true in this case for q = q x + 1. 

This completes the proof of the inequalities (48). 

From (48) we now infer at once that 4>k,i,<i (t > 0) approaches zero at 
the point x = 0. 

The inequalities (48) furnish the sufficient condition of Weierstrass 
for the absolute and uniform convergence of the series (34) in the in- 
terval < x < c. Accordingly the series (25) converges absolutely 
and uniformly in any sub-interval < e < x <[ c. 

It remains now to show that we have a solution in the system 
of functions z kl . For this purpose let us select any sub-interval 
< € < x ^ c. Since the point x = is excluded, the coefficients b'^'i 
are continuous in this interval ; and in it the system of differential 
equations (16) is satisfied by the functions z kl . For if we multiply the 
series (25) for z tJ by h'jf A we obtain the absolutely and uniformly con- 
vergent series : 

g=oo 

bXi z ij= 2 h Xi z ij, q - 

5=0 

Taking the sum of such series for all values of i anij we have : 

%—m. j=e t i=m j—e i 3=00 

(57) 2 2 #*/= 2 2 2 ft**. 

i=l j—1 t=l J=l 3=0 

3=00 j_ i=m j~=£ { 



g=oo j_ «=m J~e t - — 1 

= 2 |2 2 ^A 

3=0 L 1=1 j—1 - 1 
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where the last transformation is valid since we are dealing with abso- 
lutely convergent series. Also : 

1 g=cc 1 

(58) x Z *-'- 1= ^ x Zi ''~ x " 

3=0 

(59) ^=2°^,,,,. 

3=0 

Adding (57), (58), (59), and changing slightly the summation on the 
right, wo have : 

(60) - s,,,_ x + -* z k ,, + 22 K -' Zi -* 
X X i=i J=i 



1 r q 

= - **, /-l, + — 2*. /, + ^ I ^ **• '-!• «'+ 1 + ~ **. '• S-+ 1 
* X s=0 I- 

i—m j. 



o=0 \~ X 

2 2 J ^, • 
i=i j=i -i 



+ 



The series on the right is an absolutely and uniformly convergent 
series of continuous functions. If we replace the terms by their values 
given in (18) and (19), we can also write the series in the form : 

d ^ d "^ d 

(61) d* **•<•<> + Z dx Zk ^ +1 ~ % Tx Z "- ! «> 

3=0 3=0 

which is the series of derivatives of the terms of (25). From this it 
follows that, if we differentiate the series for z k j term by term, we shall 
obtain an absolutely and uniformly convergent series of continuous 
functions ; and therefore z k t has a continuous derivative at each point 
of our sub-interval, which is precisely (61), and this, as we have seen, 
is the same as the right side of (60). We can therefore write (60) in 
the form 

d 1 r h i=m i=ei ij 

di Zk - 1 = x Zkj - 1 + i Ha + 2 2 h ' k < iZ u> 

t=i j=i 

and now giving k and I all possible values we have precisely the system 
of equations (16). We have, then, the following result: 

In any sub-interval of c, not including the point x = 0, the functions 
z ta represented by the series (25) are continuous in x, have continuous 
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first derivatives with regard to x, and satisfy the canonical system of 
equations at every point. The terms of the series for z kl are given by the 
formulae (26), in which the functions 4>k,i, q are continuous in x through- 
out the whole of the interval b, and reduce to zero when x = 0, except 
when q = 0. The functions 4>i,i,o are constants, zero in all cases but the 
following : 

§5. 

Linear Independence of the Solutions of the Canonical 

System. 

We have shown that corresponding to each elementary divisor 
(r — r K ) C/t , there are e« solutions of (16) obtained by giving A the values 
1, 2, . . . e K ; and for the development of these solutions we have re- 
quired that | b'jfj | | log x \ e K~ x shall be integrable up tox= 0, where 
e K is determined for the root r K by the condition (23). We have, then, 
n solutions which may be written as follows : 

( when Rr h i$. R r K 
z"' = x r " 4>*l <ot Rr k — Rr K , k ± k and I <J L, 

(or k = k and I <| A, 

(62) *"' * = x* (log x) l ~ L tf * Rr k = Rr K , k dp k a.nd L < I < e k , 
Z K ' l =x-(logx) 4,^ \<l<e K . 

K A 

where the functions <$> ' are continuous in a; and 

« a, C ^ i « 

™*,j to " ( or k — k. and £ < A, 

(63) J 

<iUo=(r="X)T A<^*, 

( K =l, 2, ...i») (A. = 1, 2, . . . O. 

It is worth while to note three facts in regard to the z's, which will be 
useful later on : 

K A 

I. 2 ' does not involve log x explicitly whenever \>l; 

(64) II. limit af"' /' * = when Rr K > Rr,; 

III. limit x - r " (log a:) - e "~ a£ =0 in all cases except the 



1=0 
one, k = k and / = e K ; and then the limit is 



(«« - X) 1 ' 
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We will now show that the n solutions (62) are linearly independent. 
Suppose they were not and that there were n constants (7 Kj a not all zero 
such that : 
/c-\ X? ^C n «•* n [k = l,2,...m\ 

K=l A=l \ > ' K/ 

It will be convenient to suppose that our notation is such that : 

(66) Rr t <; R r 2 < . . . < R r m . 

Consider first those equations of (65) for which Rr k = R r 1# We 
have: 

k— m A— e K 

(67) limit ^ 2 C ^ x ~ rk *U = ° *'* = **■ 

- 1 *=* (Z=l, 2, ...**). 

Now let I = 1 in (67), and consider the limit of each term for any 
given value, within the range indicated, for k. For each term in which 
Rr K ~> R r k — Rr x the limit is zero by II. We have left, now, only 
the terms : 

(68) C^k*****.). Rr K = Rr k = R ri . 

According to I. no logarithms appear explicitly in (68), and we can 
write : 

( 69 ) C^ k *Ti = G .J^ k *tv 

Now by (63) the limit of all such terms is zero except in the one case 
k — k and A = 1, and for this term we have : 

(70) limit C kfl x- r « *£J = limit <7 M <&\ = C kA . 

So in the case of I = 1 the limit (67) turns out to be O k i when 
we evaluate the limit term by term. Now this is impossible unless 
G k<1 = 0; and so we must have, writing now k in place of k: 

(71) C„ 1 = Rr K = Rr 1 . 

Now consider in the same way the cases of (67) in which 1—2; and 
choose any one of the values of h indicated. Here again by II. the limit 
of each term is zero when i? r K > Rr k ; and we have left the terms : 

(72) C ; A aT* £ * Rr K = Rr k = R n 

\>2. 
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Now by T. no logarithms appear explicitly in the terms (72), for the 
only case in which they could occur would be for X = 1, but by (71) 
such terms do not appear. We can then write each term of (72) 



73) C K>K x- r «zti= 0,, a * r ''- r * *k A ; 

and by (63) the limit of each is zero except in the one case : 

(74) limit C k 2 x~ rk **■ } = limit <? M <ft » = Q. , . 

z=0 ' x=0 

We have then, reasoning as before, and including the previous result, 

(75) C«,i=C« i s = Rr K = Rr v 
Now the same reasoning can be applied until we have finally 

(76) C«,a = 0; Rr K = R ri , 

X> 1 

and there are left in (65) only those solutions for which Rr K > Rr^ 
Having thus disposed of all the r's in (66) such that Rr K = Rr x , we 
consider the next set of r'a whose real parts are all equal and as small as 
any other in the new set of r's, and show, in exactly the same way, that 
the corresponding constants, C Kj \ are zero. Continuing in this way, we 
finally reach the result that all the constants in (65) are zero, and that 
the assumption made as to the dependence of the n solutions leads to a 
contradiction. The n solutions are therefore linearly independent. 

§ 6 
Kkturn from the Canonical System to the Original System. 

We will now return to the original system of equations (1). Each 
solution of the canonical system (1 6) determines a solution of the original 
system (1) on account of the relation (10). We have then the following 
n solutions of (1) : 



(77) #*= 2 2 A 



i.k.l*k,l 



,K,\ 



k=l 1=1 



i = 1, 2, . 


.n 


K = 1, 2, . 


. . m 


X=l,2, . 


. e K 



Now these n solutions are linearly independent, for the z solutions are 
linearly independent, and the determinant of the constants A ikjl is not 
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Now consider the following limit : 

(78) limit af T '< (log x)~ (e '~ h) yf • x 

= limit 2 2 ^***~ r "(U>gaj)~* , «- A) aJ l \ 

a=0 fc=l fc=l 

Ou the right, if we take the limit of each term in the summation, we 
shall find that it is zero in every case except the one in which k = k and 
/ = e K (cf. (64) III.). We have, then : 

(79) limit *-*■" (logs)-*^*? * = A«,«« (g|e A X )I ' 

So we can write as a set of solutions of (1), corresponding to the ele- 
mentary divisors (r — r K ) e,L : ' 

(80) y t - =x*{\o g x)* K (^ =1;2) _J, 
where the functions i/£' K are continuous in the interval < x < c, and 



such that : 




(81) 


^ Uo- (g K _x)! J, '' t ' e < 



The constants .4,;*,/ are determined independently of the functions 
a, ^ in (1) ; and therefore all that we have said on page 346 in regard 
to certain sets of them as linearly independent solutions of the equations 
(4) in the special case of a {J = holds equally well here. 

In order to obtain solutions corresponding to (r — r K ) e «,* we assumed 
in the treatment of the canonical system that |4^| |log x\ eK ~ was inte- 
grate up to x = 0. Now, since the coefficients Vj£ t are linear functions 
with constant coefficients of the coefficients a id , it will be sufficient, in 
order to obtain solutions of (1) corresponding to (r — »"„)**, to assume 
that \a,-j\ | log xf K ~ is integrable up to x = for all values of i andy. 

Our results may be stated as follows : If (r — r^f K is an elementary 
divisor of A(r), and if we consider all the elementary divisors (r — r 4 ) e * of 
A(r) such that ftr k = Rr K , and denote by e s that exponent which is as 
great as any other exponent in this set of elementary divisors, and assume 
that | a,j\ | log x f K ~~ is integrable up to x = ; we can develop e K solutions 
of(l): 

* Cf. p. 349. 
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where i/.j' A is continuous in the neighborhood of x = cmd 

If the root refurnishes s elementary divisors : 

(r-r K )\ (r-r K ) e <+\ (r - r .)'•*— l , 

then the constants : 

-at,K,e K ) Ai tlc ^.i <eic _^ l , Ai t k+s—1, e K _|_ s _i (t = 1, 2 . . . «) 

are s linearly independent solutions of the equations (4) when r = r<. 

7%e « solutions of the differential equations (1) ite we obtain when 
k = 1, 2, . . . »» are linearly independent. 

§7. 

The Homogeneous Linear Differential Equation of the 

wiA Order. 

We shall consider homogeneous linear differential equations which can 
be written as follows: 

(82) ^ + \x +pl )d^ + \lc' + P *Jxd^ + --- 

in which pi, i>a, • • . p„ are constants, and pi, Pt, • . • p n are functions 
of the real independent variable x, continuous in the interval < x ^ b, 
and such that their absolute values are integrable up to x = ; in short, 
these p's are to have the same properties as the functions a {J in (1). 

This equation can be replaced by a system of linear differential equa- 
tions by the following substitutions : 

(83) x'- 1| = 9lM (.' = 0, 1, . . . n - 1). 
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We thus get the system of differential equations : 

(84) 



d y* - 1 „ a. 



Vi 



(i = 2, 3, . . . «). 



The characteristic determinant of (84) is : 



(85) 



r — n + 1 + ^ ix<i 

— 1 r — n + 2 











-l 





r — 


1 





— 


1 


r 



The adjoints of the first line are : 

(86) r(r — 1) ... (r — n+2) ; r(r—l) ... (r — n + S);... 

r (r — 1) ; r ; 1. 
The characteristic equation is, then : 

(87) »• (r — 1) . . . (r — n + 1) + H r (r — 1) . . . (r — n + 2) + . . . 

+ A<„-2 r (r — 1) + ,*„_! r + /*„ = 0* 

There is always one first minor of (85) which is not zero, the adjoint 
of it,, ; and therefore if r K is a multiple root of (87), the only elementary 
divisor of (85) corresponding to r K is (i »«)% where e K is the multi- 
plicity of the root r K . 

We have seen that, corresponding to the elementary divisor (r — r K ) e * , 
there are e K linearly independent solutions of (84) of the form (80). 
Here the constants A^ K>eic must satisfy the set of equations : 

(88) (r. - n + 1 + Ml ) A^ Ke< + ^ A 2>K<e< + ... + /,„ 4,,.,^ = 

-4-1, «,.« + ('■.-» + •■) 4, « 1<i( = o 

(« = 2,3, ...«). 
This system of equations has essentially only one solution, namely : 



(89) 



A,K,e K = pr t (r K — l)...(r K — n+l + i) (i = 1, 2, . . . n — 1), 

A n ,K,e K — P> 

* (87) is also called the indicial equation of (82) for the point a; = 0. 
vol. xxxviii. — 24 
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where p is a constant not zero. We may divide each solution by the 
corresponding number , and the resulting solutions may be 

\€ K A) . 

written : 

(90) j^aT- (log xy^JS^ (izllt','.'.'. ".) ' 

where E? K ^ = r K (r, - 1) . . . (r K - n + 1 + i), 
E*' x \ _ = 1. 

By means of (83) we can now return to solutions of the equation (82) 
with the following result: 

If r K is a root of the characteristic equation, and \pt\ | log x\ K ~ is 
integrable up to x = 0. where e K > e k for all k's such that Rr k = Rr K , 
e k being the multiplicity of the root r m then the equation (82) has e K 
linearly independent solutions which may be written with their first n — 1 
derivatives -. 

y ' =x" (log x) * E^ 

(91) ctv K ' K ■ „ > (A=l, 2, . . .e«) 

dx' v ° ' "-' 

where the functions E^ t are continuous in the neighborhood of x = 
and, 

(92) " ^ |a:=0 

E n-% \x=o = r K {r K — 1) . . . (r K — i + 1). 

Even for the equation of the second order this theorem does not 
merely give the results of Professor B6cher's paper above quoted, but 
goes a step farther, since in the case in which the two roots of the 
characteristic equation are equal, we require merely that : 
I b 

I (log x)\p 1 \ dx, I (log x) \p 2 1 dx 
"o 

converge, while Professor Bocher's method made it necessary for him to 
require that 

b b 

I (log xf | p 1 1 dx, I (log a;) 2 \p 2 \dx 



converge.* 

* L. c. p. 48. The function x q x in this paper is the same as our pj. 



